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The mean-field Kuramoto model for synchronization of 
phase oscillators with an asymmetric bimodal frequency dis- 
tribution is analyzed. Breaking the reflection symmetry facil- 
itates oscillator synchronization to rotating wave phases. Nu- 
merical simulations support the results based of bifurcation 
theory and high-frequency calculations. In the latter case, 
the order parameter is a linear superposition of parameters 
corresponding to rotating and counterrotating phases. 
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Collective synchronization and incoherence in large 
populations of nonlinearly coupled oscillators received a 
great attention in the recent years. Motivation for this 
can be found in the broad variety of phenomena which 
can be modeled in this framework. Indeed, synchronous 
flashing in swarms of fireflies crickets that chirp in 
unison [01, epilectic seizures in the brain Sj, electri- 

arrays 
, some 



cal synchrony among cardiac pacemaker cells 
of Josephson junctions ||^, chemical processes [5| 
models of charge density waves in quasi-one-dimensional 
metals [0, and some neural networks used to model dy- 
namic learning processes Q, all seem to be described in 
these terms. 

The mathematical model conceived first as a large col- 
lection of elementary nonlinear phase oscillators, each 
with a globally attracting limit-cycle, goes back to Win- 
free 1^ . It was later formulated as a system of nonlinearly 
coupled differential equations by Kuramoto ||l0|| , in the 
mean-field coupling case, and as a system of Langevin 
equations, (adding external white noise sources), by Sak- 
aguchi fntl. 



N 



N 



Here, 6i{t) denotes the ith oscillator phase, Ui its natu- 
ral frequency (picked up from a given distribution gito)), 



K > represents the coupling strength, and the ^^'s are 
independent identically distributed white noises. Con- 
sequently, the one-phase oscillator probability density, 
p{9,t,uj), obeys the following nonlinear Fokker-Planck 
equation, in the thermodynamic limit N ^ 00: 



dp 
dt 



de^ 39 

where D > Q comes from the noise terms in (|^), and 



v{9, t,uj) ^ uj + Kr sui{ip - 9). 
Here the complex-valued order parameter. 
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fined by 



re 



^ p{9,t,Lo)g{uj)dujd9. 



is de- 



(4) 



It is understood that (^ must be accompanied by the 
prescription of the initial value p{9,Q,tjj) = po{9,Lu), 
27r-periodic boundary conditions, and normalization 
J^"" p{9,t,Lo)d9^l. 

The fundamental phenomenon of transition from in- 
coherence [p = l/(27r), r = 0] to collective synchro- 
nization (r ^ 0) is similar to phase transitions in Sta- 
tistical Physics, and it was first analyzed rigorously by 
Strogatz and MiroUo [|l2|. They studied the linear sta- 
bility of incoherence of populations characterized by uni- 
modal frequency distributions. In fl^ , a nonlinear sta- 
bility analysis was accomplished, and bimodal frequency 
distributions [g{i-u) with two peaks] were also consid- 
ered. In the latter case, new bifurcations were discov- 
ered, showing the existence of a rich phenomenology, 
such as subcritical spontaneous stationary synchroniza- 
tion, supercritical time-periodic synchronization, bista- 
bility and hysteretic phenomena. A large amount of in- 
formation was obtained in ||l^,|l^ , adopting as models of 
uni- and bi-modal frequency distributions, g{Lu) — 6(uj), 
and g{uj) = ^[S{uj+u!q)+6{u! — ujo)]. It may be surprising 
now to realize that the asymmetric bimodal distribution. 
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g{cj) = a5{uj - ujq) + (1 - a)6{uj + ivo), (5) 



entails essentially different features with respect to the 
symmetric case, a = ^, even for a close to ^. Since, due 
to unavoidable imperfections, possibly small deviations 
from symmetry are most likely in Nature, the asymmet- 
ric case should be rather ubiquitous. The purpose of 
this paper is to illustrate the distinctive features of an 
asymmetric oscillator frequency distribution. The main 
qualitative effect of asymmetry is that no synchronized 
stationary phase is possible. Synchronized phases branch 
off from incoherence as traveling waves (TW, see below) 
and their structure becomes richer as the strength of the 
coupling increases. Asymmetry of the frequency distri- 
bution changes the stability boundaries of the incoher- 
ence (see the phase diagrams in Figs. 1 and 2), rendering 
it less stable, and, consequently, rendering the partially 
synchonized solution (whose order parameter, however, 
now depends always on time) more stable. 

The stability boundaries for the incoherent solution, 
Po{0, t^) = 1/ (27r) can be calculated by setting to zero the 
greatest of the i?e(A)'s, where p = po + ee^*^ 'r]{6, uj) (e — > 
0), and A are the eigenvalues of the linearized problem. 
They are given by 



-no Z + 1 



= 1, 



(6) 
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where we have defined x = K/ D and Z = 1 + \/ D . For 
the asymmetric bimodal distribution (|^), we find 
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.xy 



y^ + ,r^(2a_i) 



1/2 



(7) 



where y — uji^t/D. The stability regions in figs. 1, 2 are 
then determined by the condition max Re Z\_i < 1. 
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FIG. 2. Stability boundaries for the incoherent for the 
asymmetric bimodal frequency distribution, a — 0.49 (solid 
line), a = 0.3 (dotted line). 



The branch on the right of the asymptote in fig. 2 
is not completely unexpected. Indeed, its counterpart 
in the symmetric case is a parabolic profile continuing 
that in fig.l (see [^). In the latter case, however, such 
a branch is not as important as in fig. 2, since it does 
not separate different stability regions. The behavior de- 
picted in figs. 1 and 2 is confirmed by direct numerical 
simulation of the Kuramoto-Sakaguchi equation 

(^); see the evolution of the amplitude and phase of the 
order parameter in figs. 3, 4, and 5. 




FIG. 3. Time evolution of the amplitude \r{t)\ and phase 
tp{t) of the order parameter for K = 3,D = l,u!o —4, with 
a = 0.4. Note the stability of the incoherent solution. 



FIG. 1. Stability boundaries for the incoherent solution for 
the symmetric bimodal frequency distribution. Incoherence 
is linearly stable in the region to the left of the solid line. 
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FIG. 4. Time evolution of the amplitude \r{t)\ and phase 
^{t) of the order parameter for K = 4:,D = 1,ujo = 4, with 
a = 0.4. This is the region of stability of the TW solution 
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FIG. 5. Time evolution of the amplitude \r{t)\ and phase 
^{t) of the order parameter for K = 6,D = l,ujo = 4, with 
a = 0.4. Another TW solution has bifurcated from inco- 
herence resulting in a more complex behavior of the order 
parameter. 

Observe that the phase iplt) is always time-dependent, 
rather than a constant as in the case of the symmetric bi- 
modal distribution ||l^,0 . The new synchronized phases 
are described by a bifurcation analysis near the line in the 
parameter space where the incoherence loses stability: 



D 



vMr^V(!-f)(T^-f 



(8) 



This is obtained when the largest real part of the eigen- 
values is set to zero, and corresponds to Kc — 4Z?, 
Wo > D, of the symmetric case (see fig. 1). The two- 
time asymptotic analysis conducted in ||T^ may be used 



unchanged for bifurcations at the line with the asym- 
metric frequency distribution, taking into account that 
now n'^ — Lu'^ + — KcD/2, and that g{u}) is the asym- 
metric frequency distribution in (^. In fact, the sym- 
metric case possesses the reflection symmetry ujq — > — woj 
9^—9, which causes the eigenvalues to be doubly de- 
generated iQ , whereas this is not the case for the asym- 
metric frequency distribution. Then, the simple analysis 
of Ref. (which overlooked eigenvalue multiplicity, as 
pointed out in ; see also ) can be directly used for 
the asymmetric case. The result is that a branch of sta- 
ble synchronized phases bifurcates from incoherence at 
the point given by (^). Near the bifurcation line, these 
solutions have the form of TWs rotating counterclock- 



wisely |13 



1 n t-i'o{K-K,)(,t-to) 
ZTT 2n[U + i(uj + 



-m^ +C.C. (9) 

+0{K-K,), 



where c.c. means taking complex conjugate of the pre- 
ceding term, and 



R = 



' {K - Kc) ReAi 



Re7 



^'o = ImAi — Im7 



ReAi 
Re7 



(10) 



See Ref. for the explicit expressions of the parame- 
ters 7 and Ai — [dX/ dK)\K=K^- In the symmetric case, 
another solution corresponding to waves rotating clock- 
wisely has to be added to (||). This results in a stable 
standing wave solution, whose order parameter has a con- 
stant phase and an oscillatory amplitude | |l^ , p^ . 

In the high-frequency limit, loq oo, & different per- 
turbation analysis provides expressions for the evolution 
of the probability density, either near or far from bifur- 
cation points [l^ . The main result is that the frequency 
distribution decomposes in as many phases as peaks of 
the oscillator frequency distribution in such a limit. Each 
phase rotates with the frequency corresponding to its re- 
spective peak. Then, the order parameter may be written 
as a linear superposition of the order parameters of the 
different phases. For the asymmetric bimodal distribu- 
tion, the overall order parameter evolves (except by a 
constant phase shift) to 

r e*'^ = ai?+e*'^°*+*+ + (1 - a)i?_e-^'^°*+*- , (11) 

where R± and \l/± correspond to phases rotating with 
angular speeds iwo- They can be calculated with the 
stationary formule (2.1) and (1.7) of Ref. with zero 
frequency ||T^. Let a < 1/2 to be specific. We have 
the following possibilities depending on the value of the 
coupling constant: 

1. If < X < 2_D/(1 — a), the incoherent solution 
Pq = l/(27r) is stable and it is the only possible 
stationary solution. 
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2. If 2D/(1 - a) < K < 2D/a, a globally stable 
partially synchronized solution branches off inco- 
herence at K ^ 2L'/(1 - a). It has R+ = 0, 
■0 = — cjo^i and r = (1 — a) R-. Its component 
p+ = l/(27r) is incoherent, while its component p_ 
is synchronized. The overall effect is having a TW 
solution (rotating clockwisely) . 

?,. li K > 2D /a, the component p+ becomes par- 
tially synchronized too. The probability density 
then has TW components rotating clockwisely and 
anticlockwisely. Their order parameters have dif- 
ferent strengths, and i?_ > R+ if a < 1/2. 

Let us now compare the analytical results obtained in 
the high-frequency limit with those obtained by means 
of bifurcation theory. As ujq —^ oo, the parameters Ai,7, 
and in (M) become (cf. p|), 



Ai = 



2D' 



(l-a)ir2 



7 
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We can now calculate the order parameter in 
®; a-iid the previous expression: 



re'i' « (1 - a) 



K{1 



D 



-icjQt 



by using 



(12) 



Eq. ( p^ ) agrees exactly with the results of the high- 
frequency limit (|Tl|) in ||l^ : The amplitude of the order 
parameter is constant, and its phase decreases linearly in 
time. Of course, for larger values of the coupling con- 
stant another branch of oscillatory solutions (TW rotat- 
ing clockwisely) bifurcates from incoherence. Then, the 
overall probability density is richer, with an order param- 
eter whose amplitude and phase both vary with time as 
in Fig. 5. 

In conclusion, we have analyzed the mean-field 
Kuramoto-Sakaguchi model of oscillator synchronization 
with an asymmetric bimodal frequency distribution. In 
this case, reflection symmetry is broken, which results in 
stable synchronized phases that have the form of TWs 
(rotating clockwisely or anticlockwisely). These waves 
have order parameters with constant amplitude, and 
phases which depend linearly on time. As the strength 
of the coupling constant increases, such a synchronized 
phase bifurcates from incoherence. Larger values of the 
coupling strength result in a new bifurcation, which con- 
tributes to another TW. Then, both phase and amplitude 
of the order parameter become time-dependent. Numer- 
ical simulations of the model favorably agree with the 
results of bifurcation theory, and of high-frequency per- 
turbation expansions. Extensions of our analyses to the 
case of a multimodal frequency distribution [i.e., a dis- 
crete, or a continuous g{i-u) having m peaks] are worth 
considering in future works. 
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